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We present symmetry properties of the lattice vibrations of graphene nanoribbons with pure 
armchair (AGNR) and zigzag edges (ZGNR). In non-symmorphic nanoribbons the phonon modes 
at the edge of the Brillouin zone are twofold degenerate, whereas the phonon modes in symmorphic 
nanoribbons are non-degenerate. We identified the Raman-active and infrared-active modes. We 
predict and 'i{N -\- 1) Raman-active modes for A^'-ZGNRs and A-AGNRs, respectively [N is the 
number of dimers per unit cell). These modes can be used for the experimental characterization of 
graphene nanoribbons. Calculations based on density functional theory suggest that the frequency 
splitting of the LO and TO in AGNRs (corresponding to the E2g mode in graphene) exhibits 
characteristic width and family dependence. Further, all graphene nanoribbons have a Raman-active 
breathing-like mode, the frequency of which is inversely proportional to the nanoribbon width and 
thus might be used for experimental determination of the width of graphene nanoribbons. 

PACS numbers: 61.48. De 63.22.-m 63.20.D- 63.20.dk 



I. INTRODUCTION 

Graphene, a single layer of carbon atoms, attracted 
tremendous attention for both theoretical and exper- 
imental studies. Its unique properties as a two- 
dimensional crystal make it a fascinating model system 
for fundamental studies of condensed matter physics, 
while its remarkable electronic and transport properties 
and the available fabrication by lithographic processes 
are promising for future application in nanoelectronic de- 
vices, e.g., ballistic room-temperature transistors. How- 
ever, defect-free and unpatterned graphene is a zero-gap 
semiconductor and is thus unsuitable for use in novel 
electronic devices that require precise control over carrier 
type and transport. In this context, graphene nanorib- 
bons (GNRs), narrow stripes of graphene, have evolved 
from model systems for the investigation of edge effects 
in graphene and graphit^ into promising materials in 
their own right. While those nanoribbons possess simi- 
larly outstanding properties as graphene, edge effects like 
the quantum confinement of the electronic wavefunctions 
in sufficiently narrow nanoribbons can result in the open- 
ing of a band gap in the electronic structure, exhibiting 
a characteristic dependence on chirality and nanoribbon 
width^"^. Furthermore, they give rise to interesting mag- 
netic properties in nanoribbons with zizag-shaped edges 
(ZGNRs)i6 8 13^ 

Lattice vibrations have a significant effect on the trans- 
port of valence electrons in graphene-related materials 
due to a strong electron-phonon coupling^^. Therefore, 
the investigation of vibrational properties of graphene 
nanoribbons is of great interest for the physical under- 
standing of those structures. So far, there are only a few 
reports addressing this topic, the most of which focus 
on edge-localiz ed ph onons in armchair graphene nanorib- 
bons (AGNRs)li^. Kawai et al.^^ and Malola et al^ 
predict the occurance of a characteristic Raman peak at 
a frequency about 2000 cm~^ for nonhydrogenized arm- 
chair nanoribbons due to edge-related vibrations. Zhou 



et alW^ calculated the lattice vibrations at the F-point 
of various AGNRs and ZGNRs and found a Raman ac- 
tive phonon mode of width-dependent frequency, simi- 
lar to the radial breathing mode (RBM) in carbon nan- 
otubes. These insights could prove to be useful for char- 
acterization purposes. Yamada et alW^ used a molec- 
ular approach to derive the full phonon dispersions of 
small AGNRs and ZGNRs by zonefolding of the calcu- 
lated phonons of polyaromatic hydrocarbons (PAHs) of 
suitable sizes and thus were able to relate molecular vi- 
brations and phonons in crystals. 

However, to the best of our knowledge, no systematic 
work on the symmetry properties of the phonons in GNRs 
has been reported so far. The use of crystal symmetries 
is a versatile instrument in solid state physics, e.g., for 
the theoretical investigation and understanding phonon 
dispersion, electronic band structure and optical activity, 
which are ruled by symmetry and selection rules. In this 
sense, not only edge effects but also the different symme- 
try in graphene nanoribbons compared to graphene and 
carbon nanotubes (CNTs) could be responsible for differ- 
ences in electronic and vibrational properties. Moreover, 
symmetry-based selection rules for optical methods, such 
as Raman and IR spectroscopy, will help to predict and 
understand experimental results. 

In this work, we present the symmetry properties of 
graphene nanoribbons with pure armchair and zizgag 
edges. We use group theoretical methods to derive the 
full dynamical representations of AGNRs and ZGNRs, 
describing the full symmetry properties of their phonon 
spectra, and discuss Raman and infrared active modes. 
Further, we compare our results with ah initio density 
functional theory calculations of the nanoribbon phonons 
from previous worli^il. 
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II. CALCULATIONS 

The calculations of plionon frequencies of GNRs were 
performed with the computational package SIESTA'221211^ 
utilizing a density functional theory approach in the local 
approximation forrrPH. Pseudopotentials were generated 
with the TrouUier-Martins schemd^. The valence elec- 
trons were described by a double-^ basis set plus an ad- 
ditional polarizing orbital. The localization of the basis 
followed the standard split scheme and was controlled by 
the energy shift, an internal SIESTA parameter, which 
we set to a value of 50meV. Integrations in real space 
were performed on a grid with a fineness of 0.08 A, which 
can represent plane waves up to an energy of 270 Ry. 30 
fc-points along the reciprocal lattice vector were used to 
approximate integrations in reciprocal space. The space 
between periodic nanoribbon images in our case was at 
least 20 A in order to prevent interaction between them. 
We fully relaxed the atomic positions of both AGNRs and 
ZGNRs until atomic forces were less than 0.01 eV/A. The 
phonon calculations were performed by applying the fi- 
nite difference method^''. We used a supercell approach 
with a 9 X 9 X 1 supercell, the above parameters for cut- 
off energy and basis set and a grid of 3x3x1 k-points 
to calculate the phonon dispersion of graphene. All fre- 
quencies were scaled by a constant factor to match the 
experimental frequencies of the i?2g-mode in graphene, in 
order to achieve a better comparability between theory 
and experiment. 

III. RESULTS AND DISCUSSION 

A. Symmetry Properties 

In this section, we discuss the determination of the irre- 
ducible representations that correspond to the phonons of 
a given nanoribbon. As graphene nanoribbons are quasi- 
ID-crystals and possess translational periodicity in only 
one direction, the description of their symmetry proper- 
ties can be done using line groups^^M] j^^g ^^le pe- 
riodicity of the crystal, all combinations of point group 
symmetry elements, i.e., rotations, mirror planes and in- 
version symmetry, with translations that leave the crys- 
tal invariant, are elements of the crystal symmetry group 
and thus elements of the line group L. The action of a 
symmetry transformation 

^ {Q,\t + € L (1) 

on an arbitrary vector r of the nanoribbon is then a trans- 
formation of r by the point symmetry operation Qi with 
a successive rigid translation by (r -I- j/j)a, i.e., 

{Qi\T + i^i)r = Qjr + (r + i^i)a 

where r is an arbitrary integer and denotes a translation 
by a multiple of the lattice vector a and i^i e [0; 1) 



denotes a translation by a fraction of a. 

A phonon disturbs an underlying crystal by displacing 
the crystal atoms by a small amount. The representation 
of the phonons r^''"" is thus given by the direct product 
of the equivalence representation T*^', which mirrors the 
symmetry properties of the nanoribbon, and the repre- 
sentation r"*^"^ of an arbitrary polar vector (representing 
the displacement of the phonon) . It can be decomposed 
into irreducible representations, i.e., 

YPhon _ peg ^ puec 
j 

where aj is the number of each irreducible representa- 
tion in r^'''°". Only those lattice vibrations are normal 
modes which are compatible with the irreducible rep- 
resentations of r^'*"". The symmetry properties of the 
phonons with wave vector k can be extracted by decom- 
posing the equivalent representations T'^^ and T^'"^ and 
applying Eq. [2j The characters of the equivalence repre- 
sentation for the symmetry operation Ri are given by 

i 

where the delta function 5 is 1 for atoms that are shifted 
by the symmetry operation Ri from the position Vj to 
the position of an equivalent atom, and otherwise. The 
exponential function in the sum adds a phase factor for 
phonons with wave vectors k, for which Qik = k -I- K, 
where K is a reciprocal lattice vector. In case of nanorib- 
bons, K = — e, for k = - and for all other k vectors, 
including the F-point. 

1. Phonons at the T -point 

The symmetry properties of graphene nanoribbons de- 
pend on whether the number of dimers per unit cell 
N is even or odd. As visible in Figs. [T] and [2] the 
point group components Qi of the symmetry operations 
{Qi\vi+t) e L of all nanoribbons are three two-fold rota- 
tions C2,C2,C2 and three mirror operations (t^*',ct^^,(t^^. 
This includes a point inversion / = cr^^Cf. Therefore, 
the factor group, i.e., the group of symmetry operations 
Ri = {QiWi)j neglecting the translations by integer mul- 
tiples of the lattice vector a, is isomorphic to the point 
group D2h- 

For iV-AGNRs with odd N and iV-ZGNRs with even 
N, Vi can be set to zero for all symmetry operations Ri , 
see Fig. [l] (b) and Fig. [2](b). These nanoribbons belong 
to the symmorphic line group L2/mmm. Only atoms 
transformed onto equivalent atoms by the symmetry op- 
erations (QilO) contribute to the equivalence representa- 
tion and, using the character tables for the line group 
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FIG. 1: (color online) Symmetry operations {Qi\T + ly) 
for armchair nanoribbons with (a) an even number N of 
dimers per unit cell (8-AGNR) and (b) odd number N 
of dimers (7-AGNR). The operations {E\t) and (cr^^|t), 
which transform each atom in the unit cell (thick grey 
lines) onto itself or an equivalent atom in another unit 
cell, are not shown. Arrows denote twofold rotations 
about the x-,y- and z-axes, dashed lines are mirror or 
glide planes. Dimers are emphasized by red (thin grey) 
lines. 



L2/mmm (see appendix p3] Table IV ), we obtain 

N 



-peq, k=0 

^ A, odd N 



-{oA+ ® 0^ ) ® ^^TT^ioAt ® ) 



N{oA+®oA+) 



^1=0 — qAq 



) oAt © oAl 



where Ta and Tz denote representations for AGNRs and 
ZGNR, respectively. The line-group notation for the ir- 
reducible representations follows Refs. [2S] and [23] 

The atomic displacements of corresponding atoms in 
different unit cells are equal for phonons with fc = 0. In 
this case, we can express the dynamical representations 
of the phonons by means of irreducible representations 
of the point group D2h (For conversion of the line-group 
notation into the molecular notation of irreducible rep- 
resentations of point groups, see appendix [b| Table IV I. 
Eq. [2] then yields 



-pphon, k=0 
A, odd N 



N{Ag ® Bi„ ® B2g ® Bsu) 
^^^(A,®Big) 



(4) 




(O'' IX+0.5)*(C2lX+0.5) 



(CJIX) 



FIG. 2: (color online) Symmetry operations {Qi\T + 
for zigzag nanoribbons with (a) an odd number TV of 
dimers per unit cell (5-ZGNR) and (b) an even number 
N of dimers (4-ZGNR). As for Fig. [l] {E\t) and (cr^^^lr) 
are not shown. The same notation as in Fig.jljis used. 



As all components for both AGNRs and ZGNRs are 
one-dimensional, there are no systematically degenerate 
phonon modes at the F-point. 

In case of armchair nanoribbons with even N and 
zigzag nanoribbons with odd N, an additional frac- 
tional translation with v — 0.5 is needed for some 
symmetry operations Ri to maintain crystal symmetry. 
These nanoribbons thus belong to the non-symmorphic 
line group L2i/mcm (appendix |b] table |V| . In case 
of armchair nanoribbons, the equivalence representation 
changes compared to odd-iV-AGNRs due to the fact 
that the nanoribbon axis does not contain any carbon 
atoms, see Figll] (a). This results in XA'ovcnN(C2) = 
X7:t^°Ni'^^'') = 0- In ZGNRs however, "the altered 
atomic structure has no effect on r'"''''='', as the char- 
acters do not change. The equivalence representations 
are then 



N{oA+®oA, 



"pcq, k=0 "pcq, k = 

Z. odd N ^ Z. even P 



TIZ, = qA„ ® aA+ ® oA^ 

Applying Eq. [2] results in the dynamic representation 
for the F-point phonons of A^- AGNRs with even N and 
ZGNRs with odd N: 



r^r^erN = N{Ag ® Bl„ ® B2g ® Bsn) 

N 

®^iAu(S Big ® B2u ® Bsg) 



pphon, k = . 

^ Z, odd N 



N{Ag ® Blu ® B2g I 

Big ® B2u) 



B: 



3ii 



(6) 



(7) 



= N{Ag ® Bi„ ® B2g ® B^u 
® Big ® B2u) 



(5) 



Note that the equivalence representations for odd-A^- and 
even-A^-ZGNRs are equal, which implies that the symme- 
try properties of phonons at the F-point are the same for 
all ZGNRs. 
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2. Phonons with wavevector k 

The same approach as for phonons at the F-point is 
apphcable for general lattice vibrations as well, with dif- 
ferent representations for the phonon wavevectors Y'"^'^''^. 
We derived the full dynamical representations for the 
phonons of nanoribbons with armchair and zigzag edges, 
respectively, which are listed in appendix [B] 

The wavevectors k and -k are equivalent and thus rep- 
resentations for k S (0,-) and —k G (—-,0) are two- 
dimensional representation's!!^ which have to be compat- 
ible with the representations at the F-point. These com- 
patibility relations can be used to construct the repre- 
sentations of the phonons with 0<|fc|<^. This is done 
by merging those representations of the phonons at the 
F-point that are equivalent under the lowered symmetry 
for 0<|fc|<^. The displacement of the nanoribbon atoms 
due to a phonon with a wavevector k and a frequency uj 
is given by 

cj,{x, y, z, t) - A{x, y, Xo,yo, ^o, ^o) exp^^''^'-^o)~^^t-to)) 

(8) 

where xq, yo, zq are the coordinates of the origin and to 
is the initial time. Eq. [8] indicates that rotations around 
the X- and the y-axis [(Cf |t -I- i^cx), (Cf 1^ + ^Cy)l a-nd 
the reflection through the xy-plane {a^y\t + v^^y) are 
no symmetry operations for phonons with wavevectors 
0<|fc|<^. Also, those phonons obviously possess no in- 
version symmetry {I\t-\-v), i.e., the parity of the F-point 
representations vanishes when they are extended into the 
Brillouin zone. As a result of the lowered symmetry and 
the equivalence of k and -fc, the representations qAq and 
o^Q , which are equivalent for the symmetry operations 
{C^\t + vcz), {(7''^\t + i^axz) and {a^^lt + i^ayz), induce 
the two-dimensional representation ~^Ea„- In the same 
way, the representations qA^ and qA^ induce the repre- 
sentation ^\Eai- qBq and o^i^ induce the represen- 
tations '^^Esg and ~^Ebj^, respectively (refer to Sec- 
tion |A|. Thus, the dynamical representations can also 
be derived from compatibility relations to the F-point 
phonons instead of using Eq. [2j 

For symmorphic nanoribbons. the compatibility rela- 
tions can be used for phonons at the edge of the Brillouin 
zone. There, the representations are one-dimensional 
and equivalent to the eight representations at the F- 
point. Pairs of the two-dimensional representations for 
k = (0, ^), ~\Exy, transform into one dimensional rep- 
resentations of the same type, one representation of even 
parity {X^) and one with odd parity {X~), respectively. 
Thus, there is no systematic degeneracy of phonon modes 
at the edge of the Brillouin zone for symmorphic nanorib- 
bons. 

For non-symmorphic nanoribbons, pairs of represen- 
tations ~\Eao and ~\Eai {~\Ebo and ^^Eb^) induce 
two-dimensional representations ^-E^^ (--^Sq) at fc = 
^. Consequently, in non-symmorphic nanoribbons, all 
phonon modes are twofold degenerate at the edge of the 
Brillouin zone. 



The symmetry-induced degeneracy of the phonon modes 
at the Brillouin zone can be clearly seen in ah initio cal- 
culations of the full phonon dispersions of small-width 
nanoribbons presented in previous worli^. We found 
that almost all branches of the phonon dispersion of AG- 
NRs with even N and ZGNRs with odd N, i.e., the 
non-symmorphic nanoribbons, are twofold degenerate at 
fc = ^ . No systematic degeneracy was found for the sym- 
morphic nanoribbons, i.e., iV-AGNRs with odd N . Cor- 
responding results have been reported by Yamada et aP. 
These results were not understood from simple zonefold- 
ing considerations. Here we have shown that instead the 
degeneracy of phonon modes at the Brillouin zone edge 
is solely determinded by the symmetry properties of the 
nanoribbons. We refer to Fig. 9 in Ref. [2T and Figs. 3 
and 9 in Ref. ||20j for calculated phonon dispersions of 
various AGNRs and ZGNRs. 



B. Raman active modes 

In this section, we will analyse which of the nanoribbon 
phonon modes can be observed in a Raman experiment. 
Raman scattering is inelastic scattering of light, where a 
change in frequency of the scattered light occurs due to 
the creation or the annihilation of Raman active phonons 
in the material. The polarization and the intensity of 
the scattered light are linked to the incident light by the 
Raman tensor TZ; the scattering intensity is proportional 
to 

|e, •7^•e,|2 

where and e<; denote the polarization of incoming and 
scattered light, respectively. As the Raman tensor is a 
second-rank tensor, its symmetry r'^"'""'" is given by the 
direct product of the representations of a polar vector, 

pRaman ^VeC,k—Q ^ ^VeC.k — Q* 

= Ag® Big ® B2g ® Bsg (9) 

Only the phonon components contained in the represen- 
tation of the Raman tensor are Raman active. All Raman 
active phonon modes are even (subscript g) under inver- 
sion. 

The polarization of the scattered light in relation to the 
polarization of the incident light depends on the involved 
Raman active phonon, see Table |T] In the case that a 
phonon with Ag symmetry is emitted or absorbed, the 
polarization of the scattered light is parallel to that of 
the incident light. In case of Big, incoming and scattered 
light are cross-polarized in the xy-plane. Similarly, i?2g 
and B^g change the polarization of the incoming light in 
the a;z-plane and yz-plane, respectively. 

In a former work, we showed that it is possible to 
classify the F-point phonon modes of nanoribbons into 
six fundamental modes and 6(iV — 1) overtones, — 1 
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TABLE I: Allowed Raman scattering configurations. 



phonon symmetry 


allowed scattering configuration (e^, e^) 


Ag 


{z,z) 


Big 




B2g 


{x,z) 


Big 





overtones for each fundamental mode. The six funda- 
mental modes are similar to the six phonon branches 
at the F-point of graphene, thus they are of special in- 
terest. The symmetry properties of the acoustic funda- 
mental modes are equal for all AGNRs and ZGNRs, and 
none of these modes is Raman active. The longitudinal- 
acoustic phonon (0-LA), which induces an atomic dis- 
placement parallel to the nanoribbon axis, has Biu sym- 
metry. The representations of the transverse-acoustic 
(0-TA) and the out-of-plane-acoustic fundamental mode 
(0-ZA) are B^u and -6211, respectively, for all nanorib- 
bons. In contrast, the representations of the three op- 
tical fundamental modes are of even parity, and con- 
sequently, these phonons are Raman active. For AG- 
NRs, the longitudinal-optical fundamental mode (0-LO) 
has Ag symmetry and can thus be observed only for zz- 
polarization. The transverse-optical (0-TO) and out-of- 
plane-optical modes (0-ZO) have i?2g and B^g symmetry, 
respectively, and are forbidden in zz-polarizaton. We 
therefore expect that the three optical modes can be ob- 
served separately by polarization-dependent Raman mea- 
surements. 

For ZGNRs, the representations of the 0-LO and 0-TO 
are interchanged compared to AGNRs {i.e., 0-LO: i?2g, 
0-TO: Ag) and the 0-ZO has Big symmetry. The reason 
for this is the difference in structure between AGNRs and 
ZGNRs. In ZGNRs, the twofold rotation Cf around the 
X-axis of the nanoribbon and the mirror plane axy leave 
carbon atoms invariant, in constrast to AGNRs. In this 
sense, there is a simple correspondence between ZGNRs 
and armchair nanotubes, where the TO is fully symmet- 
ric and Raman-allowed in zz-configuration, and between 
AGNRs and zigzag nanotubes, where the LO phonon is 
fully symmetric. Note that in AGNRs and ZGNRs both 
LO and TO are Raman active, though in different scat- 
tering geometries, whereas in zigzag and armchair car- 
bon nanotubes the TO and LO, respectively, are Raman 
inactiv^^. For characterization of nanoribbons or de- 
termination of the graphene crystallographic direction, 
the different LO and TO selection rules in AGNRs and 
ZGNRs might be useful. 

As shown in Ref. lU there are three families of arm- 
chair nanoribbons regarding the F-point frequencies of 
the 0-LO and the 0-TO. This corresponds to the be- 
haviour of the electronic properties'"' , where AGNRs were 
found to be quasi-metallic or semiconducting, depending 
on the number of dimers in the unit cell. The magni- 
tude of the band gap can be classified into distinct fam- 




Nanoribbon width w (A) 

FIG. 3: (color online) Difference in frequency between 
the transverse-optical (0-TO) and the 
longitudinal-optical (0-LO) fundamental modes in 
armchair nanoribbons. The size of the LO-TO-splitting 
is noticeable different for AGNRs of different families. 

TABLE II: Symmetries of the Raman active overtones 
in AGNRs and ZGNRs. Due to inversion symmetry, the 
Raman active overtones of optical fundamental modes 
are all of even vibrational order. Similarly, the Raman 
active 'acoustic' overtones possess an odd vibrational 
order. 







n even 






n odd 






n-LO 


n-TO 


n-ZO 


n-LA 


n-TA 


n-ZA 


AGNR 


Ag 


Big 


Big 


B2g 


Ag 


Big 


ZGNR 


B2g 


Ag 


Big 


B2g 


Ag 


Big 



ilies. This family behaviour is also found for the LO- 
TO-splitting, see Fig. [3j The quasi-metallic nanoribbons 
of the (iV=3p+2)-family possess the largest splitting of 
the 0-LO and the 0-TO, which indicates a Kohn-anomaly 
at the F-point for these nanoribbons. The frequency 
splitting quickly decreases with increasing nanoribbon 
width w. In comparison, the nanoribbons of the (N=Sp)- 
family, which are semiconducting with intermediately 
sized band gaps, exhibit a similar but weaker drop in fre- 
quency splitting. (3p+l)-AGNRs, possessing the largest 
band gaps, exhibit only a weak dependence of the LO- 
TO-splitting on the nanoribbon width. We found an 
overall small decrease from 14 cm^^ for a 4- AGNR to 
11 cm-i for a 22-AGNR. 



1. Raman activity of overtones 

An inspection of the symmetry properties of overtones 
shows that all overtones with an even vibrational order 
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possess the same symmetry properties as their respective 
fundamental modes (at the F-point), e.g., 0-LO, 2-LO; 
4-LO, etc. in an AGNR are all of Ag symmetry. All 
these modes have in common that there is no vibrational 
"node" coinciding with the nanoribbon axis, see discus- 
sion in Ref. El] 

In the same way, the overtones with an odd vibrational 
order of a particular fundamental mode possess the same 
irreducible representation, and there is always a vibra- 
tional node on the xz-plane of the nanoribbon. A major 
point regarding the Raman activity is the inversion sym- 
metry of the phonons. The acoustic fundamental modes 
and their overtones of even order are antisymmetric with 
respect to inversion and are thus not Raman active. The 
vibrational node on the nanoribbon axis in case of odd- 
order overtones of the acoustic modes, however, induces 
a inversion symmetry. Consequenty, those phonons are 
Raman active, see Table |ll] For the optical phonons, 
overtones of even vibrational order possess even parity 
under inversion. However, the overtones with odd vi- 
brational order possess odd parity under inversion and 
consequently are not Raman active. 

A nanoribbon with an even number N of dimers per unit 
cell has ^ — 1 overtones of even vibrational order and ^ 
overtones of odd vibrational order for each fundamental 
mode. In case of a nanoribbon with an odd number of 
dimers, there are ^^^^ overtones of even and ^^-^ over- 
tones of odd vibational order. Thus, we would expect to 
find iN Raman active phonons in case of nanoribbons 
with even N and 3(A^-i- 1) Raman active phonons in case 
of nanoribbons with an odd number N of dimers per unit 
cell. The Raman active modes and their symmetries are 
summarized in Table |lll 

Regarding spectroscopic characterization we note that 
the B^g symmetry is found only for AGNRs for the 
out-of-plane-optical fundamental mode and its overtones. 
This could provide a way to experimentally distingiush 
armchair and zigzag edges in nanoribbons by cross- 
polarized Raman measurements, see Tables |T] and |lT] 
However, it will be experimentally difficult to perform 
a measurement with light polarization perpendicular to 
the nanoribbon plane. 

2. Breathing-like mode in nanoribbons 

Another promising phonon mode for characterization 
purposes is the first overtone of the transverse-acoustic 
fundamental mode (1-TA). Here, all atoms of one half of 
the nanoribbon move in-phase and in opposite direction 
than the atoms in the other half, as shown in Fig. |4] (b) . 
This results in a breathing-like expansion and compres- 
sion of the nanoribbon and shows great similarities to the 
radial breathing mode (RBM) of carbon nanotubes. The 
displacement pattern of this phonon mode is symmetric 
under inversion, thus it is Raman active. The symmetry 
of the mode is Ag (refer to Table |lT]) . From a theoretical 
point of view, its frequency should exhibit a characteristic 
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FIG. 4: (Color online) Frequency of the breathing-like 
mode (BLM) in armchair (red squares) and zigzag 
(green circles) nanoribbons of different widths w. 
Zonefolding of the linear part of the LA branch in 
graphene results in an approximation for the BLM of 
<^BLM = where a=1025.6 A- cm~^ is the sound 
velocity in graphene. ction of the BLM phonon. The 
nanoribbon is uniformly expanded and compressed. 

dependence of the nanoribbon width. Fig.|4](a) shows the 
calculated frequencies of BLMs of AGNRs and ZGNRs 
of various widths. The frequencyies of the breathing-like 
modes display a strong dependence on the nanoribbon 
width and is nearly independent of the edge type of the 
nanoribbon, i.e., the BLM has almost the same frequency 
for AGNRs and ZGNRs of equal widths. Small deviations 
are found for the smallest nanoribbons. 

We assigned wavelengths and wave vectors to the 
BLMs of various nanoribbons and mapped them onto 
the phonon dispersion of graphene. We showed in a pre- 
vious work that the T-point phonons of AGNRs can be 
unfolded to reproduce the FKM direction in graphene, 
while vibrations of ZGNRs reproduce the phonon dis- 
persion in FM direction, respectively. In this way, the 
BLMs of nanoribbons of increasing width should move 
along the LA branch of graphene towards the F-point in 
the phonon dispersion. Figure [5] shows such a mapping 
of the BLMs onto the phonon dispersion of graphene. It 
can be seen that the BLMs of AGNRs and ZGNRs exhibit 
deviations from the LA branch for smaller nanoribbons, 
but comply very well in the case of larger nanoribbons. 
The agreement of our calculated BLM frequencies with 
the linear part of the graphene LA branch enables us to 
derive a formula for estimation of the BLM frequency of 
larger nanoribbons. In the linear part of the branch, the 
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FIG. 5: (Color online) Mapping of the BLM 
frequencies of nanoribbons of various widths onto the 
phonon dispersion of graphene. As the BLM is the first 
overtone of the nanoribbon TA fundamental mode, the 
corresponding wave vector decreases with increasing 
nanoribbon width. Due to different orientation of the 
nanoribbon and the graphene lattice, the nanoribbon 
TA corresponds to the graphene LA. The BLM 
generally lies on the graphene LA branch, except for 
nanoribbons of very small widths (=larger wave 
vectors). The calculated nanoribbon frequencies have 
been rescaled by the same constant factor as the 
graphene frequencies. 



dispersion is given by 



27r 

T 



where fc is the wave length of the vibration and a = 
1025.6 A- cm ^ is the sound velocity in our calulated 
phonon dispersion of graphene. The wave length of the 
BLM is A = 2w for all nanoribbons, where w is the 
nanoribbon width. The frequency of the BLM of not 
too narrow nanoribbons can then be estimated by 



^BLM 



W 

3222A- 



(10) 



In this zonefolding model, the BLM of a nanoribbon 
with a width of 10 nm would then have a frequency of 
approximately 32 cm""'^. We expect that the agreement 
of the BLM frequencies with the LA branch of graphene 
further improves for nanoribbons of increasng width. It 



is thus likely that the Eq. 10 can be used for estimating 
the BLM frequencies of nanoribbons with realistic widths 
w > 10 nm. 



C. Infrared activity 

Another spectroscopic method for determining vibra- 
tional properties is the exitation of phonons by absorb- 
tion of infrared light. For a phonon to be infrared active. 



TABLE III: Symmetries of the infrared (IR) active 
overtones in AGNRs and ZGNRs. As all nanoribbons 
possess inversion symmetry, only phonons that are not 
Raman active are infrared active. 







n odd 






n even 






n-LO 


n-TO 


n-ZO 


n-LA 


n-TA 


n-ZA 


AGNR 










B-iu 


B2u 


ZGNR 






B2u 




Bsu 


B2u 



it has to change the permanent dipole moment in the 
solid and thus be able to couple with an electromagnetic 
field. The representation F™ of IR active phonons thus 
transforms like the representation of a polar vector r"'^'^, 
i.e.. 



i?9 



3u 



(11) 



The u indicates that the decomposition of F™ consists 
of only "odd" irreducible representations. Consequently, 
for nanoribbons, all F-point phonons that are not Ra- 
man active are IR active, i.e., the three acoustic funda- 
mental modes and the respective overtones of even vibra- 
tional order and the odd overtones of optical fundamen- 
tal modes. Thus, we would expect to find 3iV IR active 
phonons for all nanoribbons, see Table [Till 



IV. CONCLUSION 

We used group theory to study the symmetry prop- 
erties of graphene nanoribbons with pure armchair- and 
zigzag-type edges and derived the full dynamical repre- 
sentations for these nanoribbons. The different symme- 
try properties of the phonons of nanoribbons with sym- 
morphic and non-symmorphic symmetry groups mainly 
manifest themselfs at the edge of the Brillouin zone. 
Here, all phonon modes of non-symmorphic nanorib- 
bons are two- fold degenerate, whereas the phonon modes 
of symmorphic nanoribbons are always non-degenerate. 
These theoretical derivations explain results from previ- 
ous ab initio calculations^, where we calculated the full 
phonon dispersions of small-width AGNRs and ZGNRs. 
They also explain the differences in degeneracy at the 
edge of the Brillouin zone between graphene nanoribbons 
and carbon nanotubes. Using the dynamical representa- 
tions at the F-point, we predict the phonon spectra of 
graphene nanoribbons to consist of 3A^ Raman active 
and 37V infrared active vibrational modes. Some Ra- 
man active modes are promising candidates for experi- 
mental characterization of nanoribbons. Our calculations 
of phonon frequencies using density functional theory 
suggest that the frequency splitting of the longitudinal- 
optical and the transverse-optical fundamental mode 
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(corresponding to the TO and the LO of graphene at 
the F-point) exhibits a characteristic family dependence, 
with quasi-metallic nanoribbons possessing the largest 
splittings. Further, we found a Raman active breathing- 
like mode (BLM) with a strong width-dependence, sim- 
ilar to the radial breathing mode in carbon nanotubes. 
Our results might prove to be useful for future exper- 
imental investigations of the vibrational properties of 
graphene nanoribbons. 



V. ACKNOWLEDGEMENTS 

We thank C. Thomsen for his support and helpful dis- 
cussions. This work was supported in part by the Clus- 
ter of Excellence 'Unifying Concepts in Catalysis' coor- 
dinated by the TU Berlin and funded by DFG. 



Appendix A: Character tables 

Here, we show the character tables for the irreducible representations of the two symmetry groups that were relevant 
for this paper, the symmorphic L2/mmm group and the non-symmorphic L2/mcm group. Both character tables were 
adapted from previous work by Bozovic et alr^'. The corresponding matrices can be found in the same work. 



TABLE IV: Characters of the irreducible representations of the symmorphic group L2/mmm for the k-points 
k = [0; -]. For < fc < -, representations for k and —k generate two-dimensional representations ^\Ex, whereas 
the representations for fc = and fc = ^ are one-dimensional. 5q = 2cos(fcra). The first column shows the relevant 
basis functions for each fc for the purposes of this paper. The irreducible representations of the symmetry group 
L2/mmm for different k-points appear in two different notations in the second column. On the left side is the line 
group notation (LGN), on the right side is the molecular notation (MN), which is only applicable for phonons at the 
F-point. 



Basis 


LGN 


MN 


{E\r) 


(Cf|r) 


(C||r) 


(C||r) 




(o-f'lr) 


(a-ir) 


{a^y\r) 






A, 


1 


1 


1 


1 


1 


1 


1 


1 


z 




Biu 


1 


-1 


-1 


1 


-1 


1 


1 


-1 






Big 


1 


-1 


-1 


1 


1 


-1 


-1 


1 






Au 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


X 


oA+ 


Bsu 


1 


1 


-1 


-1 


-1 


-1 


1 


1 




o^r 


B2g 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


y 




B2u 


1 


-1 


1 


-1 


-1 


1 


-1 


1 






Bsg 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


z 






So{r) 








5o{r) 
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— fc 77^ 

k^Bo 




So{r) 
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Mr) 


Mr) 





X 






So{r) 








Mr) 





Mr) 


So{r) 





y 


— fc 771 
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Mr) 





So{r) 
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i-iy 


(-1)- 


z 
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(-1)- 
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-i-iy 


(-1)- 








i-iy 


(-1)- 


(-1)- 


{-ly 
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-i-iy 


-(-!)- 


X 






i-iy 


(-1)- 


-(-!)- 


-i-iy 
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-(-ly 
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y 






i-iy 


-(-!)- 


(-1)- 
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TABLE V: Characters of the irreducible representations of the non-symmorphic group L2i/mcm for the k-points 
k — (0; As for the symmorphic counterpart (table IV), the representations for < fc < ^ are two-dimensional, so 



are the representations for k — ^, which are generated by 'lE^o, \Ea-^ and \Ebq, \Eb^- = 2cos(fcra) and 
Si = 2cos(fc(r + i)a). 



Basis 


LGN 


MN 




(Cf|r) (C||r+1) 


(C||r + i) 


(i|-r+ ^) 


(cr!'^|T-+ i) 
V 1 2 ' 


\ \ J 






U 


Ag 


1 


1 


1 


1 


1 


1 


1 


1 


z 




Blu 


1 


-1 


-1 


1 


-1 


1 


1 


-1 




oB+ 


Big 


1 


-1 


-1 


1 


1 


-1 


-1 


1 




oBq 


Au 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


X 


oA+ 


Bsu 


1 


1 


-1 


-1 


-1 


-1 


1 


1 




oA^ 


B2g 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


y 




B2u 


1 


-1 


1 


-1 


-1 


1 


-1 


1 
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1 


-1 
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•50 (r) 
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y 


k^Bi 
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2{-ir 

















2{-\y 
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Appendix B: Dynamical representations of the phonons in armchair- and zigzag-edged graphene nanoribbons 

Here, we report the full dynamical representations of the phonons in ideal armchair- and zigzag-edged nanoribbons, 
as follows by Eq. [S) 

rr=°dd =^(oAj © © o^r © 0^^) © ^^^(oS+ ® o^r) © ^^^(oSq- e 0^0+) 

© Y.{2N{-\Ea, © T\Ea,) © (iV - \)-\Eb,, © + \r\EB,} 

k 

® N{kA+ ® ® kA+ ® nA^f ) ® ^^-4^{^B+ ® kB^) ® ^^^{^B+ ® kB^) (B1) 



r^even =^(o^J © O^Q^ © O^F © 0^^) © ylo^Q © qB^ © © oSf) 



k 



-fc 



IEbo © -IEbJ © N{2^E2'^ ® 



Bl^ 



(B2) 



rrcvcn =^(o^J © 0^,7 © o^r © © 0^0^ © oB+) (BNY,{'^'IEa„ ® ^-^Ea, ® "t^So © "t^sJ 



® iV(^A+ ® ^ ® ^ A+ ® ® ^B+ ® 
r^Tdd ^N{oA+ ® ® o^r © © 0^0^ © oS+) ®Ny^{2-lEA, ® 2-^^^, 



fc 



-fc JTl 

fc^Bo 



(B3) 



(B4) 
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